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ABSTRACT 

Fundamental information about the Universe is encoded in anisotropies of the Cosmic Mi- 
crowave Background (CMB) radiation. To make full use of this information, an experiment must 
image the entire sky with the angular resolution, sensitivity, and spectral coverage necessary to 
reach the limits set by cosmic variance on angular scales > 10'. Recent progress in detector tech- 
nology allows this to be achieved by a properly designed space mission that fits well within the 
scope of NASA's Medium-class Explorer program. An essential component of the mission design 
is an observing strategy that minimizes systematic effects due to instrumental offset drifts. The 
detector advances make possible a "spin chopping" approach has significant technical and scientific 
advantages over the strategy used by COBE, which reconstructed an image of the sky via inversion 
of a large matrix of differential measurements. The advantages include increased angular resolu- 
tion, increased sensitivity, and simplicity of instrumentation and spacecraft operations. For the 
parameters typical of experiments like the Primordial Structures Investigation (PSI) and the Far 
InfraRed Explorer (FIRE), we show that the spin-chopping strategy produces images of the sky 
and power spectra of CMB anisotropies that contain no significant systematic artifacts. 
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1 INTRODUCTION 



The study of the Cosmic Microwave Background (CMB) entered an exciting new phase with 
the detection of large scale anisotropics by the COBE-DMR (Smoot et al. 1992). Degree-scale 
fluctuations have now been detected by several ground- and balloon-based experiments (Cheng 
et al. 1995; de Bernardis et al. 1994; Ganga et al. 1994; Gundersen et al. 1995; Hancock et al. 
1994; Netterfield et al. 1995; Piccirillo et al. 1996; Ruhl et al. 1995; Tanaka et al. 1996). Full 
realization of the enormous scientific potential of the CMB, however, will require a high-sensitivity, 
high-angular-resolution study of a large fraction of the sky from a single, stable platform. Due 
to dramatic advances in detector sensitivity and stability over the last decade, such a study is 
now feasible, and several groups have devoted substantial effort to developing plans for a second- 
generation space experiment capable of extracting the maximum available cosmological information 
from the CMB. 

An important element in the design of a CMB experiment is the scan strategy, which in addi- 
tion to specifying how different parts of the sky are observed must provide extremely good control of 
potential sources of systematic errors. The aforementioned advances in detector technology coupled 
with recent advances in mission design allow a direct-imaging strategy that has been adopted by 
two proposed space missions: the Primordial Structures Investigation (PSI) and the Far-InfraRed 
Explorer (FIRE). This strategy has many advantages over the differential chopping scheme used 
by COBE-DMR. The purpose of this paper is to demonstrate, through analytical calculations and 
simulations, that this direct-imaging strategy can produce images of the sky containing no signif- 
icant systematic artifacts, especially those due to instrumental fluctuations with a 1/ f spectrum. 
We begin by introducing the scan strategy and its advantages. 

PSI and FIRE each propose to use a spinning spacecraft located near the Earth-Sun L 2 point, 
with the spin axis pointed at the Sun and the boresight of a single telescope at right angles to the 
spin axis. The individual beams of a focal plane array of detectors sweep out great circles through 
the ecliptic poles. The spin axis is precessed either in small steps (6.'67 nine times daily for PSI) 
or continuously (FIRE) to maintain Sun-pointing. The output of the detectors is sampled rapidly 
(> 3.5 times per beam) around each circle. An image of the whole sky is built up steadily over 
6 months as the Earth and the spacecraft orbit the Sun. Since each great circle is scanned many 
times in succession, instrumental drifts or offsets on timescales long compared to the spin period 
are easily removed. The north and south ecliptic pole regions tie any given great circle to all others. 

In contrast, the differential chopping scheme used by COBE-DMR (Boggess et al. 1992) mea- 
sures differences in power between two horn antennas or telescopes pointed at widely-separated 
positions on the sky. The horns spin around a bisecting axis, which itself is steered around the 
sky at a slower rate. After essentially the entire sky has been scanned the image is constructed 
by inverting a large matrix of temperature differences (Janssen & Gulkis 1992). The differential 
chopping strategy reduces the effects of gain drifts in detectors as well as time- varying instrumen- 
tal or atmospheric offsets, but at the price of decreased sensitivity and increased complexity of 
instrumentation and observing strategy. Historically this approach has been required when stabil- 
ity could be achieved in no other way. A variant of this strategy has been chosen for the Microwave 
Anisotropy Probe (MAP: Wright, Hinshaw, & Bennett 1995), another proposed second-generation 
CMB mission. 

The spin-chopping strategy used by PSI and FIRE has many advantages, either by itself or 
as an facilitating factor in the mission design. A complete discussion of these advantages would 
require full specifications of the spacecraft and mission designs, which are beyond the scope of the 
paper, but we introduce sufficient detail here to motivate the discussion that follows. The main 
advantages are: 
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(1) Simplicity of the hardware. Only a single telescope is required, and the complications of 
internal switches or correlation receivers are avoided. 

(2) Superior angular resolution. A single-telescope design allows the largest possible aperture 
within the size limitations of the launch vehicle. Detailed theoretical investigation of CMB 
anisotropies shows unequivocally that even modest increases in angular resolution in the critical 
sub-degree region can have dramatic scientific benefits. 

(3) Increased sensitivity. The low instrument temperatures that can be achieved passively with 
the Sun-pointed spin axis (next item) lead directly or indirectly to colder and more sensitive 
detectors. Moreover, an unswitched system has a factor of y/2 greater sensitivity than a 
switched system that differences pairs of pixels on the sky (Janssen & Gulkis, 1992). 

(4) Extremely benign and stable thermal conditions. With the spin axis pointed at the Sun (and 
< 5° from the Earth) the drivers for spin-synchronous temperature changes are essentially 
eliminated. A solar panel at one end of the spacecraft normal to the spin axis shades the rest 
of the spacecraft from the Sun and Earth at all times. Passive cooling of the instrument end 
of the spacecraft to below 60 K is readily achieved. 

(5) A direct image of the sky is formed incrementally as it is observed. One does not have to 
wait until a large fraction of the sky (for COBE-DMR almost the entire sky) is observed to 
form an image. This significantly reduces the science risk of an early mission failure, and 
allows immediate investigation of mission performance and potential systematic errors at the 
microkelvin level. 

(6) It gives a direct measure of 1// noise in the instrument, since each great circle is "closed." 
The differential strategy is vulnerable to drifts in the switched system offset. While this was 
relatively easy to handle in COBE-DMR, the problem becomes more severe at higher angular 
resolutions and sensitivities where the beams are not swapped precisely on the sky in the 
measurement of each pixel pair. Also, thermal gradients are more difficult to control with two 
large telescopes than with two small horns. 

(7) Data from different parts of the sky are independent. There is no aliasing of signal from 
high foreground regions such as the Galactic plane into foreground-free regions as there may 
be in differential chopping schemes (see also Lineweaver et al. 1994). Similarly, the effects of 
transients, either in the instrument or on the sky (e.g., variable or moving sources such as 
planets) are localized, and the large number of scans over a given great circle provides many 
cross-checks on a pixel-by-pixel basis. 

(8) Simplicity of spacecraft operations. The almost stationary spin axis simplifies attitude control 
and ground communications. 

These advantages would be lost if medium-timescale instrumental fluctuations that could not be 
"chopped out" by the spacecraft spin introduced serious systematic errors. Such fluctuations might 
be caused by, for example, detector gain variations or temperature changes in the optical system 
with a 1// spectrum. Naively one might expect that the detectors would have to have essentially 
no such fluctuations on timescales shorter than the spin period. In this paper we derive analytical 
expressions for the effect of 1// noise on sky images obtained with the scan strategy described. 
We also verify the analytical results with simulations of sky images at very high resolution, using 
detector characteristics measured in the laboratory for both transistor amplifiers (for PSI) and 
bolometers (for FIRE). We show that the naive expectation that detector fluctuations must be 
on timescales longer than the spin period is far too restrictive, and that a spin-chopping mission 
in the benign L 2 environment with recently developed and characterized transistor amplifiers or 
bolometers can realize all of the above advantages. 
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2 THEORETICAL ANALYSIS 



The spacecraft spin acts as an ideal switch or chop to suppress slowly varying drifts and 
offsets; however, the effects of instrumental drifts on timescales less than the spin period can- 
not be suppressed completely by the spin. In an instrument based on either indium phosphide 
high-electron-mobility transistor amplifiers (PSI) or AC-biased bolometers operated from a tem- 
perature regulated heat sink (FIRE), these drifts are measured to have a 1// spectrum dominated 
by transistor gain fluctuations. 

In this section we derive analytic expressions for the effect of such drifts on an image of the 
sky on a single scan circle (combining different scan circles is discussed in the next section). The 
main effect turns out to be an increase in the rms noise in the image that depends on the spacecraft 
spin rate and on the frequency at which the 1/f and the "white" noise have the same amplitude, 
known as the 1/f knee frequency. 

Assume a time-varying signal T(t) due to stationary noise that is characterized by the noise 
power spectrum S(f), where / is frequency. The signal is measured as a uniform series of contiguous 
averages over integration time At, 

Ti = / T(t) dt, (1) 

Jt + (i-l)At 

where to is some arbitrary reference time. We wish to determine the variance 

a\j-k) = l -((T ] -T k f) (2) 

in terms of the time-lag index j — k, the integration time At, and the noise power spectrum S(f). 
Angle brackets denote ensemble averages. Defining n = j — k and using the stationarity property 
of the noise, this equation may be rewritten as 

a\n) = - (T 3 T k ) . (3) 

Using Eqn. (1) and 

/•OO 

(T(t)T(t - r)> = / S(f) cos (2tt/t) df, (4) 

J 

we find 

a2{n) = 7^f S{f) ^/r^ Sin2(mr/At) df - (5) 

We have measured the detected-power spectrum S(f) of the InP HEMT amplifiers and bolome- 
ters that will be used by PSI (Seiffert et al. 1996; Gaier et al. 1996b) and FIRE, respectively. For 
both, S is well-represented by 

S(f) = a + b/f (6) 

for which Eqn. (5) becomes 



a 2 (n] 



2At 



1 + ^At <j)(n) 



(7) 



where 

(f){n) = (n - l) 2 ln(rc-l) - 2n 2 Inn + (n + l) 2 ln(n + 1) 

(8) 

— >21nn + 3 as n — > oo. 

(Convergence is rapid. For n > 10 the error is < 10~ 3 .) The logarithmic behavior of (j), essential in 
what follows, is illustrated in Figure 1. Note that the noise characteristics of many detectors cannot 
be represented by Eq. (6). For example, the detectors used on IRAS exhibited strong memory 
effects, which introduced large offsets whenever a bright object was crossed. Such detectors would 
be quite unsuitable for spin-chopped observations of the CMB. 
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Fig. 1 4>(n) from Eqn. (8). The dotted line shows the approximation (f> = 2 Inn + 3, which 
gives an error < 10~ 3 for n > 10. 

The variance cJq for a total power radiometer with pure white noise (i.e., b = in Eqn. (6)) is 
well known to be 

where T s is the system noise temperature, B is the radiometer bandwidth, and At is the averaging 
time for a given sample as defined in Eqn. (1). Eqn. (7) may then be written as 

<T 2 (n) = Cq [1 + At / knee <f>(n)] (10) 

where /k nee = b/a is the frequency at which the white noise and 1// noise contributions are equal. 
The factor 

F= [1 + At / knee 0(n)] 1/2 (11) 
thus gives the increased uncertainty over the case of perfect white noise. 
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Now consider a scan circle divided into m uniform and contiguous integrations, and assume that 
the circle is scanned N times and that the N measurements at each pixel are averaged. Specifically, 
for two pixels 1 < p < q < m around the circle, these averages are 



N 



p+(i—l)m ! 



i=l 
1 N 

and T q = — ^Tq+a-i 



(12) 



)m • 



i=l 



Let us find the variance 



a 2 (n = q-p,m,N) = U(T q -T p ) 2 ) 



(13) 



Substituting Eqn. (12) into (13), expanding, and making use of Eqn. (3), we obtain 
1 1 7V_1 

cr 2 (n, m, N) = ^°" 2 ( n ) - ~ k ) {2o- 2 (/cm) - a 2 [km + n] - a 2 [km - n] j. 

Using Eqn. (10), we arrive at 



(14) 



fc=i 



a 2 (n, m,N) = ^r [l + At / knee $ (n, m, N) 



where 



N-l 



$(n, m, iV) = 0(n) — — (iV — fe) |20(A;m) — 4>[km + n] — 0[/cm — n] |. 



fc=i 



The factor 



F(n,m,N)= 1 + At / knee $(n, m, N) 



1/2 



(15) 



(16) 



(17) 



gives the increase in image uncertainty in the scan-averaged image over the white-noise case. Fig- 
ure 2 shows <&(n,m, N) for typical parameters. The case N = 1 consists of only one scan with no 
scan averaging, so that the factor F increases simply as the two-point variance. As ./V increases $ 
becomes symmetric around the midpoint, reflecting the fact that pixels beyond the midpoint move 
closer to the first pixel of the next scan, shortening the effective time between pixel pairs. PSI will 
average 1600 scans before redirecting the spacecraft spin axis, so that we are practically concerned 
only with the asymptotic form for large N . FIRE will have N S> 1, but much smaller than for PSI. 
Note that for N > 10, $>(n,m, N) depends hardly at all on N, and a 2 oc iV -1 . 

The uncertainty is a maximum between diametrically opposed pixels if N is large. From the 
above we may write 



F max (m, N) = F(m/2, m,N) = 1 + At f knee 3> max (m, N) 



1/2 



(18) 



N-l 

$maxK N) = <t>(m/2) --^2(N-k) faikm) - <f>[(k + 1/2) m)} - cf>[(k - 1/2) m] }. (19) 



k=i 
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Pixel Separation n 



Fig. 2 $(n, m, N) from Eqn. (16) as a function of pixel separation n for a scan circle consisting 
of m = 2000 pixels. Seven curves are plotted, for (top to bottom) N = 1, 3, 10, 30, 100, 300, and 
1600. The curves for TV = 300 and N = 1600 are almost indistinguishable. For PSI, m = 1700 
(corresponding to a 12.'6 beam) and iV = 1600. For FIRE, m = 2700 (corresponding to an 8' beam) 
and N = 18. For large m the curves are insensitive to m. 

For large N the sum may be approximated as an integral and the term in {• • •} is approximately 
4>"(x) ~ x~ 2 . The sum is thus only logarithmically dependent on m. Numerical evaluation gives 
0.871 for N > 50 almost independent of m. Hence we can approximate 

$max(m) = hm $ max (m, TV) ~ <f>{m/2) - 0.871. (20) 

N^oo 

Further, if m > 10, we can simplify this to 

$ max (m) ~ 21nm + 0.743. (21) 

The maximum excess noise factor for a typical PSI or FIRE scan circle can therefore be approxi- 
mated closely as 

F max ~ [1 + At / kn ee (2 In m + 0.743)] 1/2 . (22) 
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F max is plotted in Figure 3 as a function of /knee f° r the PSI and FIRE spin rates. Note that the 
sensitivity of a differential chopping instrument corresponds to F max = y/2. 




Fig. 3 F max from Eqn. (22) versus /knee for the PSI and FIRE spin rates of 10 and 0.1 rpm, 
respectively. The pixel size for the PSI curve (lower) is 12.'6, while that for the FIRE curve is 8'. 
The symbols mark knee frequencies measured in the laboratory for stabilized InP HEMT amplifiers 
for PSI (10 Hz; Seiffert et al. 1996; Gaier et al. 1996b) and bolometers for FIRE (0.03 Hz). 

Based on laboratory measurements, /knee — 10 Hz for the stabilized HEMTs used by PSI 
(Seiffert et al. 1996; Gaier et al. 1996b), and /knee — 0.03 Hz for the bolometers used by FIRE. 
These knee frequencies correspond to increases in the noise over the S(f) = constant (pure white 
noise) case of only about 24% for PSI and 6% for FIRE. This is significantly less even than the 
minimum increase of noise in a switching system of \/2. 
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3 SIMULATIONS 



The foregoing calculation gives the variance between pixels on a given scan circle, which does 
not depend on the mean value. The variance between pixels on different scan circles, however, does 
depend on the mean values of the two circles. Since we are interested only in CMB anisotropies, the 
mean level itself is unimportant. On the other hand, since the mean value of a 1// noise component 
formally diverges as the averaging time increases, the mean levels on different scan circles will in 
general be different, and if not removed would produce stripes parallel to the scan direction. 

Assume for the moment that the offsets between scan circles have been determined and re- 
moved. The variance between pixels on the circle is unchanged. The variance between random 
pixels on different circles depends on how accurately the offset was removed and on the relative 
phases of the dominant 1/f fluctuations along the two circles. If we assume that the offset has 
been subtracted with an error small compared to cx max = ^(-Fmax/iV) 1 / 2 , then a^ ax should be an 
upper bound on the variance calculated between random pixels on different circles. 




Fig. 4 All-sky PSI simulation at 90 GHz, including both white and 1/f noise. The FIRE sim- 
ulation at a similar frequency is indistinguishable from this one. All images are Aitoff projections 
in ecliptic coordinates, with 12.6' pixels, to match roughly the resolution of the printer. (For a 
standard printer at 600 dots per inch, this gives 1.25 dots per pixel, and the image is therefore black 
and white, not grayscale.) a) Assumed sky, using a standard CDM model to generate CMB 
fluctuations together with extrapolations of Galactic data taken at other frequencies. The stretch 
is ±300 /iK. b) Same as before, but including instrument noise (white noise + 1// noise as mea- 
sured in the laboratory at 90 GHz) and the dipole. The stretch is ±3000 to accommodate the 
dipole. c) Observed sky, following fitting and removal of the dipole, plus DC-offset and gradient 
terms in each great circle. The stretch is ±300 fiK. d) Difference between lower left and upper left 
images stretched to ±30 fiK, showing the high SNR achieved. Residual striping is seen along lines of 
constant longitude at the level expected from the calculations in § 2. Stripes parallel to the sides of 
the paper are caused by uneveness in the printer, not the experiment! Note: the figure above has 
been drastically reduced in resolution in order to be a manageable 100 k postscript file. Full resolu- 



tion and colour versions of this figure can be found at [http : //astrophysics . jpl . nasa.gov/PSI 
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To complete our demonstration that the spin-chopping strategy leads to no significant sys- 
tematic errors due to medium-timescale instrumental fluctuations, we now describe simulations of 
full-sky images that: 1) demonstrate that offsets between scan circles can be removed; 2) verify the 
conjecture that cr^ iax is a good estimator of the variance in the full-sky images; and 3) show that 
there are no significant artifacts of any kind in the images. 

We simulated a CMB signal at < V resolution assuming a standard Cold Dark Matter (CDM) 
model with h = 0.5, = 1 and £Ib = 0.05. To this cosmic signal we added Galactic emission 
extrapolated from IRAS and DIRBE (for the dust) and ground based measurements (for free-free 
and synchrotron), as well as the CMB dipole. 

A time series of signal data was generated by "scanning" the simulated sky according to the 
PSI and FIRE prescriptions, to which was added a simulated time series of noise. This time stream 
was binned into the pixel on the sky corresponding to the spacecraft attitude. The noise time 
stream was generated in runs of length corresponding to ~ 2 x 10 6 pixels. Since each great circle 
contains ~ 6 x 10 3 pixels, this allowed us to wrap 60 (30) times on each scan circle for PSI (FIRE), 
and cover 6 adjacent scan circles with correlated noise. For FIRE we simulated the actual number 
of rotations per great circle, while for PSI we stopped for practicality at N = 60. As can be seen in 
Fig. 2, we are approaching the asymptotic regime at this N, as varying the number of rotations in 
the simulations confirmed. By choosing N = 60 we were able to obtain coherence in the simulated 
time stream over several great circles in the sky to maximize our sensitivity to striping. Lower 
resolution simulations with even more correlated great circles showed no extra cross-scan striping. 

INPUT OUTPUT DIFFERENCE 




Fig. 5a PSI simulation of a small patch, enlarged for clarity. Same as Figure 4, for a 14° x 14° 
equatorial region (where sensitivity and the possible effects of 1/ f noise are worst). The pixel size 
is 3.'6 (i.e. 3.5x oversampling) . The saturation level has been set at ±200 /iK for this plot. The 
bottom row is smoothed to the 12.'6 beam size. The scan direction is vertical. Note: Full 6 



panel (8k) version of this figure available at http://astrophysics.jpl.nasa.gov/PSI/ 
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For both PSI and FIRE, the power spectrum of the simulated noise was derived from real 
amplifier data. These data supported the assumption of a uniform random phase distribution for 
the simulated data, allowing us to generate the time streams in the simulation from the known 
power spectra. The noise added into the time stream was scaled relative to the signal map to the 
level expected after foreground removal from 2 years of observation. This level was determined 
in separate simulations (not described here) of foreground removal that used the raw instrumental 
noise levels (see Gaier et al. 1996a for 90 GHz performance) and beam sizes expected at each 
frequency (similar but preliminary studies are described in Brandt et al. 1995). 

For PSI the data stream came from 44 GHz tone-stabilized HEMTs, sampled at 500 Hz for 500 s. 
We derived a power spectrum that was well-fitted by white noise plus 1/ / noise with /knee = 10 Hz. 
Similar measurements for FIRE bolometers gave a spectrum well-fitted by white noise plus 1// 
noise with /knee = 0.03 Hz. 

The N rotations were averaged together to produce one great circle of data. Offsets between 
scan circles were removed simply by subtracting the mean on each circle. (For FIRE a gradient 
was also subtracted, to remove the small residual asymmetry shown in Figure 2 for N = 30. In 
practice this has essentially no effect on sky power, since the sky must "close" around the circle.) 
The result is shown in Figure 4, which looks identical at 600 dots per inch for either the PSI 90 GHz 
channel or the FIRE 100 GHz channel. Note that time domain drift removal in our configuration 
is easier than it was for COBE-DMK (e.g. Bennett et al. 1994), since we reference to the poles in 
every scan. 

INPUT OUTPUT DIFFERENCE 




IN. SMOOTHED OUT. SMOOTHED DIFF. SMOOTHED 



Fig. 5b FIRE simulation of a small patch, enlarged for clarity. Same as before, with a 
10° x 10° equatorial region. The pixel size here is 2.'25x2.'25. The bottom row is smoothed to 
the 8' beam size. The scan direction is vertical. Note: full 6 panel (8 k) version of this figure 



available at http://astrophysics.jpl.nasa.gov/PSI/ 



Figure 5a shows a 14° x 14° patch of sky for the PSI experiment enlarged for clarity to the 
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point that the printer can show a true gray scale. The beam size is 12.'6 and the pixel size is 3.'6, 
set by the sampling rate. The top row is shown at full pixel resolution; the bottom row is smoothed 
to the beam size. Left to right the panels show the assumed and observed skies, and the difference. 
The signal-to-noise ratio (SNR) is 1.1 even on the pixel scale. There are small residual "stripes" 
along scan lines (vertical here), exactly as expected from the theoretical calculations. An additional 
destriping step was performed for the PSI small patch simulations, where the mean of each scan 
across the patch was subtracted. 

Figure 5b shows the same thing for the parameters of the FIRE experiment at 220 GHz, the 
highest frequency channel at which CMB fluctuations are likely to be usefully measured. The 
beam size is 8' and the pixels are 2.'25, again set by the sampling rate. The patch plotted here 
is 10° x 10°, and the panels represent the same things as for the PSI figure, except that here the 
lower three panels have now been smoothed by 8'. The SNR achieved by PSI and FIRE will allow 
super-resolution, i.e., information can be extracted on angular scales substantially smaller than the 
beam size, and in practice one would make images that over-sample the beam by significantly larger 
factors than we have used here. 

To quantify in a different way the effect of the residual striping along scan lines we have 
extracted the power spectrum of fluctuations from our simulated skies. Simulated full-sky images 
at the highest achievable resolution were observed with the PSI or FIRE strategy, including the 
effects of noise. From these we calculated the Cg spectrum, where Ct is the usual squared amplitude 
of the average spherical harmonic at multipole I (see e.g., White, Scott & Silk 1994). Each input 
Cg is one realization of the underlying average sky, and hence the power spectra show the effects of 
cosmic and sample variance at each multipole. We assumed for convenience that foregrounds had 
been removed over 80% of the sky with attendant increase in the noise level (see above) . The other 
20% of the sky (the Galactic plane) was ignored. In practice the fraction of the sky contaminated 
with foregrounds may be somewhat higher. The only effect of this would be a corresponding 
increase in sample variance. 

Figure 6 shows the spectrum of multipole moments derived from our input, output, and dif- 
ference images. Figure 6a is for a PSI simulation at 90 GHz, with the 12.'6 beamsize. Figure 6b 
shows the same thing for FIRE, with a beam-size of 8'. 

Each plot has two panels: the upper panel has £(£ + l)Cg on the vertical axis, so that the plot 
is power spectrum per logarithmic interval in £; the lower panel has Cg on the vertical axis, this 
being the natural way to indicate noise, since pure white noise is then a horizontal line. The C/s 
extracted from the input image are shown shifted vertically by two orders of magnitude to separate 
them from the output image Cgs. Note that both the input and output power spectra are the 
underlying theoretical spectra multiplied by the experimental window function, which is why they 
drop off faster at high £ than the more familiar theoretical curves. Also note that, as expected, the 
output power spectrum is the sum of the input and the noise spectra. 

Four points are worth emphasizing: 1) no significant artifacts are apparent in the output sky 
power spectrum; 2) the difference image spectrum (i.e., the noise estimator) is nearly white noise 
over all angular scales, with some extra power at large scales where the SNR is largest; 3) the 
variance over the sky is cr^^, as expected; and 4) we expect to be able to determine the power 
spectrum efficiently and without smoothing out to £ ~ 900 (for PSI) and £ ~ 1400 (for FIRE). 

These imaging simulations reinforce our analytic calculations. Not only are there no significant 
artifacts in the images, there is none in the power spectra either. The increase of the RMS noise over 
the pure white noise case is exactly as predicted by Eq. (22). Moreover, the simulations demonstrate 
that it is straightforward to analyze the time-stream data and produce a two-dimensional image of 
the sky. Many refinements of the method are possible, particularly concerning removal of offsets 
between scan circles, but even the simple method we used here works extremely well. 
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Multipole / 



Fig. 6a PSI power spectra. The three curves plotted are from the input, output, and difference 
images. We took one realization of a standard CDM sky and observed it with the PSI spin-chopping 
scheme. The top curve is the power spectrum of the (noiseless) input data, shifted upward to avoid 
overlap. Note that the theoretical power spectrum has been convolved with the experimental beam. 
The jaggedness of this curve is just a reflection of the statistical nature of the model. The middle 
curve is the power spectrum fitted to the output image, equivalent to the lower left panel of Figure 4. 
No smoothing has been applied in i, i.e., every individual £ mode is calculated. The bottom curve 
is the power spectrum of the difference image (i.e., output minus input), equivalent to the lower 
right panel of Figure 4. Here the beamsize was taken to be 12.'6, with 3.'6 pixels, and a noise level 
of 24 /iK for a 12.'6 resolution element. Note that the noise is close to white on all angular scales, 
and that it only dominates for I > 900. Top: the familiar way of plotting power spectra, with 
l{l + l)Ci on the vertical axis. Bottom: The same plot with Ci on the vertical axis, so that 
pure white noise appears as a horizontal line. 

4 CONCLUSIONS 

PSI and FIRE image the sky in the most direct way possible, scanning it pixel by pixel. 
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10 100 1000 

Multipole I 

Fig. 6b FIRE power spectra. The same as Figure 6a, but for the parameters of the FIRE 
experiment. Here the beamsize was taken to be 8', with 2.'25 pixels, and a noise level of 25 /jK for 
an 8' resolution element. Note again that the noise is close to white on all angular scales, and that 
it only dominates for £ > 1400. Top: the familiar way of plotting power spectra, with £{£ + l)Ce 
on the vertical axis. Bottom: The same plot with Ce on the vertical axis, so that pure white 
noise appears as a horizontal line. 

This strategy, made possible by a combination of technological advances in detectors, appropriate 
mission design, and favorable conditions at L2, has important technical and scientific advantages 
over that used by COBE-DMR, in which differences are measured between widely-spaced parts of 
the sky. The technical advantages include simplicity of the hardware and spacecraft operations, 
and benign and extremely stable thermal conditions with the Sun-pointed spin axis. The scientific 
advantages include high angular resolution, high sensitivity, and the independence of data from 
different parts of the sky. This last factor not only eliminates possible aliasing of signal from high 
foreground regions into low foreground regions and localizes the effects of transients, but also allows 
data to be analyzed as they are taken. The first slice of the sky can be seen at the first downlink 
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and analyzed immediately. Mission performance and systematic errors can be investigated almost 
immediately to microkelvin levels. 

Laboratory measurements of the transistor amplifiers and bolometers used by PSI and FIRE 
respectively show that their noise characteristics can be modelled extremely well as a combination 
of white and 1/f noise. There are no long-term memory effects, such as those familiar from IRAS, 
that would result in persistent offsets after bright objects are crossed. 

We have shown both analytically and with extensive simulations of our experimental procedure 
that 1/f fluctuations on timescales that cannot be completely removed by the spacecraft spin are 
transformed by that spin into a small increase in near-white noise that introduces no significant 
systematic errors. Moreover, we have shown that the naive expectation that /knee must be much 
less than the spin frequency is far too restrictive. PSI and FIRE can thus realize all of the potential 
advantages, both scientific and technical, of the spin-chopping strategy. 
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